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3.  STATEMENT  OF  PROBLEM 


Solitons  in  stimulated  Raman  scattering  are  a  coherent  transient 
phenomenon,  in  which  short  pulses  of  coherent  (pump)  radiation 
are  created  in  an  envelope  of  longer  wave  length  (Stokes) 
radiation . 

In  a  Raman  active  medium  with  homogeneous  broadening  solitons  can 
be  created  from  a  large  class  of  initial  conditions,  including 
those  appropriate  for  spontaneous  Raman  scattering  and  for  Raman 
amplifiers  with  weak  Stokes  probe.  The  resulting  pulse  of  pump 
radiation  will  initially  have  a  width  comparable  to  the  coherence 
decay  time. 

The  problem  addressed  in  the  work  reported  here  is  to  determine 
the  detailed  dynamics  of  soliton  propagation  in  media  with 
homogeneous  and  inhomogeneous  broadening.  Of  particular  interest 
are  the  rates  of  soliton  narrowing  and  decay  and  the  frequency 
and  temporal  character istics  of  the  developing  optical  pulse. 
Soliton  decay  occurs  in  broadened  media,  if  the  two  optical  beams 
are  not  exactly  in  resonance  with  the  Raman  transition,  a 
situation  which  to  some  extent  frequently  occurs  in  practice.  It 
is  hence  important  to  determine  the  limitations  on  soliton 
stability  imposed  by  this  effect. 

Both  numerical  and  analytical  methods  have  been  employed. 


4.  SUMMARY  OF  RESULTS 


A.  RESULTS  FOR  NARROW  SOLITONS  (HYPER  TRANSIENT  REGIME) 


If  the  temporal  soliton  width  is  much  smaller  than  the  coherence 
decay  time,  broadening  can  be  treated  as  a  perturbation. 
Analytical  results  have  been  obtained  for  this  case  by  using  the 
method  of  asymptotic  perturbation  theory,  based  on  the  inverse 
scattering  transform  (1ST)  and  also  a  method  based  on  constants 
of  motion,  which  is  more  directly  related  to  the  physics  of  the 
problem.  Both  approaches  give  the  same  results. 

The  following  results  are  obtained  for  homogeneously  broadend 
media : 

For  exact  Raman  resonance  the  width  decreases  with  gain  according 
to  an  inverse  square  root  law,  and  becomes  independent  of  the 
initial  width  in  the  large  gain  limit.  For  off  resonant  Raman 
scattering  the  soliton  amplitude,  defined  as  the  maximum 
relative  pump  intensity,  decreases  exponentially  with  increasing 
gain.  Since  a  relation  exists  between  amplitude  and  frequency 
detuning,  this  means  that  the  average  detuning  increases  in  the 
process  of  propagation.  The  decay  coefficient  is  proportional  to 
the  square  of  the  initial  frequency  mismatch.  The  width 
decreases  more  slowly  than  in  the  resonant  case,  reaches  a 
minimum  at  an  amplitude  of  0.5  and  increases  beyond  that  point. 

Numerical  studies  confirm  these  results.  The  agreement  is 
excellent  for  exact  resonance,  and  very  good  for  the  off 
resonant  case.  In  the  latter  case  at  large  gains  the  rate  of 
soliton  decay  is  found  to  be  somewhat  smaller  than  predicted, 
while  the  rate  of  soliton  narrowing  is  somewhat  larger.  The 
reason  appears  to  be  an  almost  linear  frequency  chirp,  which 
develops  between  the  leading  edge  (lower  frequency)  and  trailing 
edge  (higher  frequency)  of  the  soliton.  This  is  a  second  order 
effect  which  cannot  be  modelled  in  the  approximation  employed. 

The  effect  is  small  however,  as  long  as  the  detuning  is  small 
compared  to  the  inverse  soliton  width,  even  if  it  should  be 
larger  than  the  Raman  line  width. 

For  inhomogeneous ly  broadend  media  analytical  results  have  been 
obtained  for  the  resonant  case,  and  numerical  results  for  both 
resonant  and  off  resonant  case: 

In  the  resonant  case  tr.e  soliton  width  decreases  very  weakly  with 
gain  like  the  inverse  of  the  third  root.  This  is  due  to  the  fact 
that  for  inhomogeneous  broadening  the  frequency  spectrum  is  more 
sharply  localized.  "er.ce  the  temporal  autocorrelation  function 
has  vanishing  slope  at  the  origin,  ct.od  narrowing  occurs  only  to 

A 


second  order  in  the  width.  The  predictions  agree  excellently 
with  numerical  results.- 

In  the  off  resonant  case  two  different  cases  occur.  If  the 
frequency  detuning  is  small  compared  to  the  line  width,  weak 
soliton  decay  and  narrowing  occurs.  If  the  detuning  exceeds  the 
line  width,  however,  complicated  envelope  modulations  of  the 
optical  pulses  occur.  Soliton  amplitude  and  width  go  through 
cycles  of  decrease  and  increase,  with  an  overall  tendency  for 
decay  and  pulse  broadening.  An  analytical  theory  for  this  type 
of  behavior  has  not  yet  been  developed. 


3.  RESULTS  EOR  3F0AD  SOLITONS  (TRANSIENT  REGIME) 


In  the  transient  regime  where  the  soliton  width  is  comparable  to 
the  coherence  time  analytical  methods  based  on  perturbation 
theory  are  no  longer  adequate.  Numerical  studies  have  been 
performed,  which  show  that  the  main  features  found  in  the 
hyper  transient  case  also  occur  in  this  regime.  In  particular 
pulse  narrowing  and  decay  occurs  both  in  the  homogeneously  and 
inhomogeneous !y  broadened  case. 

The  rate  of  decay  is  found  to  be  much  less  initially  than 
expected  from  the  first  order  theory.  As  a  result  a  pulse 
narrowing  to  a  width  of  about  20%  of  the  coherence  time  is 
possible  even  in  the  off  resonant  case.  A  maximal  frequency 
offset  of  about  5%  of  the  Raman  line  width  can  be  tolerated  in 
this  case.  These  results  are  encouraging  for  soliton  experiments 
in  media  with  broad  inhomogeneous  lines  like  methane  or  optical 
fibers . 

A  detailed  discussion  of  the  results  under  sections  A.  and  B. 
above  are  given  in  the  Appendices  A.  and  3.,  which  are  self 
contained.  Appendix  A.  contains  a  mathematical  analysis  of 
perturbation  theory  for  the  homogeneously  broadened  case  based  on 
soliton  theory  (1ST).  Appendix  3.  gives  a  more  physical 
discussion  of  homogeneous  and  inhomogeneous  broadening,  which 
emphasizes  the  relavant  physical  processes,  and  employs 
differential  conservation  laws  for  photon  energy  and  momentum. 


C.  SATURATING  SOLITONS  FOR  HOMOGENEOUS  BROADENING 
( HYRERTRANSIENT  REGIME) 

For  sufficiently  high  field  intensities  medium  saturation  :  ur 
Faman  soiitons  exist  also  in  this  case,  and  the  total  soliton 
width  can  be  much  larger  than  the  pulse  rise  time,  giving  rise 


to  a  top  hat  type  temporal  profile.  Soliton  theory  is  quite 
complicated  for  this  system,  and  we  have  studied  the  resonant 
case  in  the  absence  of  Stark  shift  effects.  Detailed  results  are 
discussed  in  Appendix  C. 

If  medium  saturation  is  strong,  the  rate  of  soliton  narrowing 
per  unit  propagation  distance  is  independent  of  the  soliton 
width,  and  the  latter  decreases  linearly  with  distance.  This  is 
much  stronger  than  the  inverse  square  root  law  found  for  the 
unsaturated  case,  and  makes  this  type  of  Raman  soliton  very 
attractive  for  pulse  narrowing  techniques.  Effects  of  frequency 
detuning  and  Stark  shift  remain  to  be  studied. 
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PROPAGATION  OF  NON  RESONANT  TEMPORAL  RAMAN  SOLITONS  IN  THE 
PRESENCE  OF  COHERENCE  DECAY:  A  PERTURBATION  THEORETICAL  ANALYSIS 
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INTRODUCTION 


The  occurrence  of  solitons  in  a  Raman  medium  has  been 
demonstrated  experimentally  [1].  We  present  here  the  results  of 
an  analytical  study  of  the  evolution  of  these  solitons  when  they 
propagate  in  the  Raman  medium  m  the  presence  of  decay  of 
coherence  and  detuning  from  the  Raman  .  ine  center.  Decay  of 
coherence  (which  arises  due  to  collision  of  molecules  in  a 
gaseous  medium  i s  inc.uded  as  a  perturbation  to  the  exactly 
solvable  equations  of  transient  stimulated  Raman  scattering  in 
whi'h  this  effec-  is  ignored.  The  perturbation  appears  in  the 
equation  for  tne  ‘ime  evo.iuon  of  the  medium  variable  .  On  the 
•  ther  hand  the  phys  .  a',  pr  o.em  :s  -o  find  the  spatial  evolution 
of  tne  optica,  fields  wr.i  are  given  as  functions  of  time  at 
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The  process  of  ‘ransier.*  s-  imu  1  a  ted  Raman  scattering  is  described 
-y  the  following  equations  for  the  slowly  varying  amplitudes  for 
*  he  pump  and  the  St. -.ires  fields  A  and  A  j  .  and  the  off  diagonal 
matrix  element  for  ‘ne  polar  ::a  - i  ;r.  of  ‘he  medium,  X  [1]. 
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the  optical  pulse.  Here  and  in  the  following  we  denote  partial 
differentiation  with  respect  to  these  variables  by  a  subscript. 
The  first  term  in  equation  (3!  describes  the  decay  of  coherence 
of  the  polarization  of  the  medium  with  a  rate  e  This  term  is 

assumed  small  compared  to  the  rate  of  change  of  polarization  of 
the  medium  ar.d  is  treated  as  a  perturbation  to  the  exactly 
solvable  equations.  Equations  '  1  )  to  (3;  for  the  case  €■  -0  have 
been  solved  by  the  1ST  [2,4].  The  solutions  to  equations  ;i!  to 
(3)  which  we  shall  study  nere  are  tr.e  travelling  waves  ir.  one 
soli  ton  form.  They  are  given  by 
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The  spatial  dependence  of  q  is  given  by  a  zero  order  terir. 
corresponding  to  the  solvable  case  !  =0)  and  an  additional  term 

for  €  >0  : 
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For  the  additional  term  we  find 
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(6)  togetner  with  relations  (9)  t 
angular  variables.  For  comparison  with  earlier  work  it  is 
convenient  to  rewrite  the  result  in  tr.e  form. 
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Equations  (18)  and  (19)  can  also  be  obtained  fror.  me  f_rst  t 
constants  of  motion  for  this  problem.  The  constants  arc  me: 
particular  values  for  the  solution  (4)  to  t  6  >  are  given  by 
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The  spatial  deriva*  :ves  of  both  coefficients  ;  r.  the  presence 
the  perturbation  can  r.e  -a.  r^.ated  f  nm  equation  14’  and  g:v 
equations  equivalent  *o  IS'  and  19;  above.  This  approach  * 
be  extended  to  give  an  equation  for  me  temporal  sol  item  pcs: 
(parameter  t*  above',  whim  is  equivalent  to  equation  2 C ;  . 
temporal  dependence  of  t r.e  sol:-or,  phase  parameter  b  above 
however  cannot  be  determined  in  this  way 
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expressed  in  terms  of  *  he  parameter  and  the  relative 
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parameter  &£.  From  equation  (20)  we  obtain: 
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NUMERICAL  RESULTS 

In  figures  1  to  3  we  compare  the  results  of  a  numerical 
solution  of  equations  (1)  tc  (3)  to  the  perturbation  analysis 
above.  The  exact  one  soliton  form  of  the  optical  fields  was  used 
as  initial  condition  at  X  =C ,  while  the  medium  polarization  X 
was  set  equal  to  zero  at  r =0 .  The  temporal  position  of  the 
fields  was  chosen  such  as  to  render  the  corresponding  error 
smaller  than  one  part  in  thousand.  The  parameters  $  and  bug 
were  determined  for  the  numerical  solution  for  the  optical  fields 
fr jm  the  maximal  relative  pump  intensity  and  the  area  under  the 
temporal  pump  intensity  curve.  The  constant  A <  was  determined 
for  off  resonant  soli  tons  both  from  the  pump  intensity  and  width 
by  equation  ( 2  3  j  and  from  the  temporal  derivative  of  the  relative 
phase  between  pump  and  S-o’ces  amplitude  (see  equations  (5)  and 
i  6  ;  ;  . 

Three  different  ases  a-e  shown.  The  parameters  chosen  were  tOp  =1 
for  a*,  cases.  €  =C . I  and  Wj  =C.C  in  case  1,  £.  =0.1  and  w* 

=0.5  in  case  ?,  and  €  =0.4  and  =0.5  in  case  3.  Cases  1 

ar.d  2  satisfy  me  condition  mat  the  damping  term  oe  sma  1 1 
ompared  to  r  he  driving  term 

«  h,  v-: 

The  maxima,  propaga*  urn.  distance  was  =20.  Figures  la  ar.d  lb 
snow  me  analytical  results  for  me  square  ;t  and  for  me 

exponential  function,  equation.  '24'  ,  of  me  soliton  position, 
wmcf.  are  straight  lines.  The  dots  indicate  v.e  numerics, 
values,  which  fa.i  a. most  exactly  onto  me  ♦  n.eoret  i .  a  .  1  1  ties  . 

Figures  2a  and  2b  snow  the  same  parameters  for  case  2  wim 
detuning  A)j  =0.5.  The  agreement  is  excel  let.*  for  the  so  lit  or. 
position  and  very  good  for  the  solitor.  widtr..  Larger 
discrepancies  occur  .:.  so  1  i tor.  width  for  *  •  1 0  .  In  this  region 
the  rate  of  increase  in  sci it on  width  is  sma .isr  than  predicted 
The  reason  for  this  is  mainly  the  behavior  of  the  sol  i  *  or.  phase  . 
According  to  equations  ( 5 i  and  (S'  the  rela'ive  pnase  of  pump  and 
Stokes  field  shows  linear  t « -,e  dependence  in  *  h.e  asympto*  ir 
temporal  region  away  from  me  soliton  center  corresponding  m  a 
constant  frequency.  The  relative  phase  for  me  numerical 
solution  on  the  other  hand  shows  an  a. most  linear  frequency  chirp 
with  local  frequency  at  soliton  center  c.ose  m  the  ar.aiy*  :ra. 


result.  The  chirp  increases  with  increasing  propagation 
distance.  In  figures  2c  and  2d  we  show  the  soliton  amplitude 
(maximal  relative  pump  intensity  ^  )  and  the  constant  .  The 

latter  is  determined  both  from  soliton  width  and  amplitude  (solid 
line)  by  equation  (25)  and  from  the  relative  phase  (dotted  line) 
by  equations  ! 5 )  and  (6).  Both  methods  begin  to  disagree  at 
larger  propagation  distances. 

Figures  3a  to  3d  show  the  same  results  for  case  3,  in  which 
condition  (27)  for  the  validity  of  perturbation  theory  is  no 
ionger  satisfied.  The  numerical  results  do  however  show  the  same 
qualitative  features  as  in  case  2.  The  quantitative  amount  of 
soliton  decay  and  broadening  is  smaller  than  predicted  from 
perturbation  theory.  Figure  3d  shows  that  the  parameter  as 
determined  from  the  soliton  width  disagrees  strongly  with  the 
value  obtained  from  the  phase,  the  latter  showing  a  much 
stronger  local  detuning  (smaller  ^  )  than  predicted. 
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FIGURE  CAPTIONS 


Figure  la:  The  square  of  the  inverse  soliton  width  is  shown 

as  a  function  of  distance  for  zero  detuning.  The  initial  solito 
width  is  equal  to  1  and  the  Raman  line  width  is  equal  to  0.1.  Th 
straight  solid  line  is  the  analytical  result,  and  circles  are 
the  numerical  results. 


Figure  lb:  An  exponential  function  of  the  temporal  soliton 

center  is  shown  as  a  function  of  distance.  Parameters  are  as  in 
Figure  la.  The  straight  solid  line  is  the  analytical  result,  an 
circles  are  the  numerical  results. 


Figure  2a:  The  square  of  the  inverse  soliton  width  is  shown 

as  a  function  of  distance  for  a  detuning  of  0.5.  The  initial 
soliton  width  is  equal  to  1  and  the  Raman  line  width  is  equal  to 
0.1.  The  solid  line  is  the  analytical  result,  and  circles  are 
the  numerical  results. 


Figure  2b:  An  exponential  function  of  the  temporal  soliton 

center  is  shown  as  a  function  of  distance.  Parameters  are  as  in 
Figure  2a.  The  solid  line  is  the  analytical  result,  and  circles 
are  the  numerical  results. 


Figure  2c :  The  maximal  relative  pump  intensity  is  shown  as  a 

function  of  distance.  Parameters  are  as  in  Figure  2a.  The  solid 
line  is  the  analytical  result,  and  circles  are  the  numerical 
results . 


Figure  2d:  Numerical  results  for  the  constant  as  a 

function  of  distance  are  shown.  Parameters  are  as  in  Figure  2a 
The  solid  line  is  the  result  obtained  from  the  pumo  intensitv 
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Figure  3b:  An  exponential  function  of  the  temporal  soliton 

center  is  shown  as  a  function  of  distance.  Parameters  are  as  in 
Figure  3a.  The  solid  line  is  the  analytical  result,  and  circles 
are  the  numerical  results. 


Figure  3c :  The  maximal  relative  pump  intensity  is  shown  as  a 

function  of  distance.  Parameters  are  as  in  Figure  3a.  The  solid 
line  is  the  analytical  result,  and  circles  are  the  numerical 
resul ts . 
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function  of  distance  are  shown.  Parameters  are  as  in  Figure  3a. 
The  solid  line  is  the  result  obtained  from  the  pump  intensity 
profile,  and  circles  indicate  the  result  obtained  from  the  local 
frequency  profile. 
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ABSTRACT 


We  study  the  propagat  i  or.  cf  Raman  S'. ::  r.edia  w :  - : 
homogeneous  and  inhomogeneous  media  Ana  ly  -  :  :  a  1  me- n  ids  ism  g 
d  i  f  f  erer.- i  a  I  "inspr-a*  ion  laws  for  p:..*  -  on  energy  and  -  :ne:.’  im 
currer.-s  give  good  predirt  i  ins  for  -he  ra-es  if  sol:-  or.  r.arr  '.wing 
and  decay  in  the  hyper  t  rar.s  ier.t  regime  i;r  r»s  .  1  -  s  are  -.-er 
and  extended  into  -he  transient  regime  through  n.umer  i  rai  studies 
Inhomogeneous  broadening  is  found  -  o  lead  -  he  -  -  ®r  so'.;-  ;, 
stability  than  homogeneous  croaoeriog,  and  1  :  m  i  *  s  r.  free  ue-  -  y 
detuning  from  exac-  Parian  resonance  are  der  •.•red. 


I 


So  .  :  tons  i  r.  st  :m  ;  .area  r'ar.ar.  sra*,*er:r.g  S?S  '  wfr*  firs- 
observed  experimer.ta.  .y  :  r.  1  ?  ;  :  ;  .  St  ;-.lared  Par.  a' 

s  r  a  *  *■  » r  -  -  J  :r.v-  .  v«**-  a”  .eas"  *  w  op*  ;:a.  rpar.s  .  'r.e  -.<*•  i-  '  . 

t  r.  e  pur.p  a  r.  i  a  s  d  rear  _a.  lec  S  t  o*es  bear  a  *  a  f  r  e  ~  e  r  — 

lower  -har.  *r.e  r  ;r.:  .  7  e  ;  .  -  *=•  -- e  :  -  *  ie«.  .  ,  -  :  .a  . 

a  1 me  s -  equa  1  *  c  *  r.e  *  r  a  r.  ^  ~  r  f  r  e  : u e r.  r  v  if  a  ?  a rr. a r.  a  ,:*■;■•  <=• 

medium.  ‘rr'  .jr  wh :  r.  *  p-  paga'e  .  7~e  be*  a  :  1  «•  : 

dyr.an  :  rs  r  f  s  ■  -.  *  *  e  r  ;  j  : ^  r«-  -.  s  -  -  *  re  *  :  - e  a  •  a  1  e  ~  f  a  r  i  a  *  .  *  r. 

*  r.  e  c*  .  a  .  f  i  >  .  .  - .  .  *  *  -  a  a  .  e  is  .  .  n  c  o  rr.  p  a  r  e  ;  *  *  ■ .  e 

conerer.  :e  *.:«  c  *  *  -  e  -  a  *  •=•  d  y  a  *  a  *  e  *  e  ~  .  rr  e  and  -re  red  ; 

i  s  i  r.  t  r.6  ground  s  *  i  *  *—  '  e  *  e e  rear.  is  c»  .  w *.* s  ar.r  .if",  a  r  o 


ia-  ur.ex.c  it  eo  s*a*e 


is  in  i  .10  t*  0 ’•  i  n.  cl  s  i  ”  c*  *  r’ps  p  d  n  *  t>  d  .  w<-tv  s  an  „  i  f  *  d*.  i 

-he  pump  oear.  i  s  ever*  .a.  .y  rep  .e-«  d  If  -re  field  -  ime  3'  a  1  e  is 
comparable  -  o  -  r.e  e  •  . -e  -  rar.s  i  er  *  regime  ,  *  .e 

opposite  process  ,-e.- r  -  a  i  h  .  e  .  As  -r.e  red  i  n  reccr.es 
excited,  pr.otor.s  or.  ,<e  *  ra:.s :  er  red  f  r  on  -he  S*okes  rear,  o  -•  * 
the  pur.p,  reducing  -r.e  ar.cur.-  -f  medium  excitation  and  never  s :  r. 
the  process  of  pur.p  d->r  1  e*i-‘-  y  ■>  g  ]  . 

Ir.  the  1  i  r.  i  -  wr.ere  -he  field  time  scale  is  much  snorter 
than  the  coherence  -  ir.e  hyper-  rar.s  ien-  regime'  phc-or.s  can  be 
transferred  free-'.'  i  r.  b c  -  r.  directions.  ? a r. a r.  so  1  i  tor.s  1  ,  d  are 
a  coherer.-  -rar.sier*  :.  nenonenu.  in  which  -'r.e  oppos  i  -e 
of  pur.p  ieple  - 1  on  and  ar.pl  if  i  at  ion  occur  in  a  complet 
balanced  way,  leaving  -,.e  reliur.  ir.  the  .  1 1  ial  unexcited  =.-ate 

after  the  pur.p  pulse  has  passed.  The  result  is  a  stable 
local  iced  non  linear  wave  .  f  •-  x  i  ta  -  ion  c:  both  medium  and 
optical  fields.  which  car.  be  observed  exrer  imer.tal  Iv  as  a  pulse 
of  pump  radiati..:.  *  ."filing  thru-gh  the  medium  in  an  envelope  o 
St  Dices  radiation. 

This  phenomenon  shews  very  olo  e  ana  lories  -o  the  nh.encnenc 
of  self  induced  transparency  TIT  "  <:  ,  7  I r.  fact  an  exac* 

formal  ;c  -  res- oncer.-.-.  «=•  can  be  es  t  ab  1  i  shed  ’  letweer  the  equations 
describing  -r.e  -wo  pnenom.ena.  The  two  levels  of  -he  atomic 
system  in  TIT  or  respond  to  the  -wo  levels  of  -he  photon  ays -err. 
defined  by  -ne  -wo  beams  ir.  SP.S  ,  while  Me  elec-ric  field  in  SI 
'corresponds  to  the  coherent  polarization  ir.  S?S  .  A  1st.  -he  - jles 
of  spatial  and  temporal  variables  are  interchanged . 


Similar  pn.er.or.i-  .-t  nave  been  investigated  t  r.eore  t  i  c  a  1  1  y  for 
coherent  two  photon  propagation  ]o]  and  for  three  level  systems, 
where  all  transition-;  *  re  nearly  resonant  with  an  optical 
f  r  equen.cy  [  9  ]  . 


Soli  tons  in  S?S  have  beer,  studied  theorem 


•a  r  i  ous 


mcdel  systems  '5,10,11].  Strictly  speaking  they  exist  only  in  th 
hyper  t  rans  i  en  t  limit,  wh.ere  coherence  decay  -an  be  r.eglec-ed, 
and  the  equations  for  S?.S  are  completely  integrable.  Their 
perfect  stability  implies  also,  that  they  cannot  be  generated 


>T>  x> 


:T".er  - r.an  rrom  optica.  r:e.is.  wr.ose  tempora.  profile  c^ose.y 
T.atches  the  required  forrr.  .  These  are  of  rather  special  fore, 
which  may  be  difficult  to  generate  ar.d  precisely  time  for  the  t 
beams  involved.  and  i  t  seems  -rat  mainly  for  -his  reason  Pamar. 
soli  tons  have  attracted  little  a t tertian  before. 


C n  -  n e  o t  r. e r  nano  in  the  c r e s e r c e  of  c v - . e r e n c e  c e c a v  t  r e 
equations  are  no  .  or.ger  e:*ia<  *  1  y  i :  -eyrar  1  e  .  So  1  i  ton  like 
exc  i  tat  ions  exist  m  the  transient  and  hyper t rar.s lent  regime 
wr.  i  h.  closely  c "  r  respe-n  i  -  -cl  yis  ar.d  in  mam.v  cases  converge 
to  so.  .tins  in  -  he  limit  of  .  a  r  ge  gam  '1,12.13.14].  T.ne  i  r 
proper  -  les  cat.  oe  ar.a.yne  ■;  :  .  .decs*  t  c  i  n  terms  of  the 

integrable  sys-em  vanishing  line  width  ,  and  for  these  reasons 
we  shall  use  t.ne  -err  .  i- on  ‘  a  1  -  hough  -hey  are  not  soliton 
i  r.  -he  strict  sense  of  -he  -err.. 


These  soiuti  ms  ;ar.  he  generated  form  a  large  class  of 
initial  renditions  '1.14  i.  ,  1  h  ,  which  includes  conditi  ..r.s 
appropriate  for  Raman  amplifiers  with  small  Stokes  seed.  The 
important  ingredier.-  r.ere  is  a  rapid  phase  shift  in  the  ir.'ecte 
Stokes  beam,  which,  -riggers  -r.e  "rar.sier.t  phase  wave 
character  is- ic  for  -r.e  fu.iy  developed  sc  1  :  -  or.  .  In  fa;-  it  has 
beer,  found  ey.perimer.ta.  ly  .  14]  ar.  i  ir.  an  independent  the -.re  -  ica 

study  ]17]  -hat  soiitons  car.  even  re  rer.erafed  from  Dhase 
*  ' 

ma 

.r.  a  d  d  i  t :  o 
; le ; .  soi 
comp i  a-e 
four  wave 


the  temporal  width,  of  -  r.e  nru*=*r 

.  n  Cr035^r.Q  p  ".**  G  7;  rJ  (j*  <-d  7  .  O’.  1  .  1  ■  ^  d  T.  ^  tJ 

ar.d  S-ok.es  beam  are  no-  e:-:a*-.y 
‘  r  ar.s  .  *  .  ;  n  .  *  he  amp  f  *  1 

increasing  propagation  dis-ar.  :e 

It  is  importar.  -  - 
quantitative  unders-a-.d  i  ng  :  f  -  : 
fundamental  imperta:  «  - f  -  -  . 
respect  to  possible  appli:a-ij-. 
]23]  -ha-  soli -ary  w.»--e*,  f  -he 
- c  accelerate  -r.arged  par-  i  *  les 
•eloci-ies.  Addi-ior.ai  ar.r  .  .  ,  • 


in  rn is  paper  we 
on  the  propagation  of 
i  nhomoger.eous  i  y  fcroade 


WTVT 


differential  conservation  laws  for  energy  and  momentum.  w- : 
give  formal  expression  to  the  balance  of  physical  pro  >ss<>s 
the  system.  In  the  presence  of  broadening  -he  *  c  r  resp o  :  no- 
quantities  are  no  longer  conserved,  but  satisfy  a  system  :f 
differential  equations.  If  the  time  dependence  of  tr.e  p- : 
fields  is  assumed  to  oe  of  one  soliton  form  wi-.n  posit;;:, 
dependent  parameters,  a  simple  closed  set  of  equ 
parameters  results,  which  is  exact  to  .owest  or 
broaden i ng  .  This  method  is  closely  related  - c  and  r.as  beer, 
employed  in  the  context  of  asymp-o-ic  perturbation  -heory  case 
on  the  inverse  scattering  -nans form  IS”  [2'.. 22, 21]. 


In  the  present  context  this  approach  is  useful  from  sever 
points  of  view.  Firs- ly  accurate  analytical  results  are  or  -  a : 
in  a  simple  and  transparent  way,  without  having  -o  emplry 
powerful  and  mathema*  i  :al  ly  complex  roc*s  of  1ST.  Se'-cr.d  1  y  *  r. 
conservation  laws  and  -he  derived  dif feren-  ial  equations  nave 
direct  physical  In-^rpr^ta*  :  on  i  r.  -eras  of  -he  perfer-  -r 
imperfect  balance  of  opposite  pnysical  processes  .  abscrct  ;  a 


emission,  gain  and  loss  which  governs  the  situation  a1 


Perfect  balance  is  maintained  in  the  acser/e  of  ~ o n e r e n  ■  e  deca 
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Here  and  Af  are  the  Stokes  and  pump  electric  fields 

Q  is t  the  off  diagonal  matrix  el ement  for  the  molecular  Raman 
transition,  and  may  be  considered  as  an  amplitude  of  coherent 


:ne  mec: 


excitation  of  the  medium  'ep  ical  phonon  1  .  x  and  ^  are  time 
like  and  space  like  coordinates,  which,  are  related  to  time  t 
and  propagation  distance  z  in  the  laboratory  frame  by  T  =  t 
z  c  and  'i*  =  z  Partial  differentiation  with  respect  to  these 
coordinates  is  indicated  ry  the  corresponding  subscript.  The 
first  term  in  '2.1  describes  collisions!  coherence  decay  with 
decay  time  1  jp  wnere  f  is  the  angu.ar  Raman  line  width  HWHM 
Lorer.t  z  ian  .  m  radians  per  unit  of  time.  Suitable  units  have 
oeer.  chosen  to  render  all  coupling  constants  equal  to  unity. 


Certain  ef fee's  are  neglected  in  these  equations.  In 
addition  to  higher  order  Stores  generation  and  four  wave  mixing 
these  are  effects  of  dynamic  Stark  shift ,  medium  polarization 
and  medium  saturation  population  of  the  upper  level :  [ 1C  .  18 , 29' 
T.oeir  va  .  idi'y  is  r.er.e  .  ir.'ec  to  sufficiently  low  in  tens  it  i 

I*  is  convenient  to  introduce  a  Biccr  vector  S  for  tr.e 
o p - i c  a  1  fields  by 

« 

S^  =  S.  -  :  S,  =  2  A_,  Aq  !  2 

-  A3  A3  -  A(,  A,. 

« 

3.  -  A  _  A  -  A _  A  _ 


.  r.  terms  ~  f  s 
reformulated  as 


S3  S,  =  S^  *  53  S3 


the  equations 


z  .~  e  a  . 


The  last  equation  in  '2.5  is  a  differential  conservatio 
aw  expressing  conservat  i  or.  of  photon  number.  In  the  absence  of 
o.oerenr  e  decay  i  r -0  :  equations  .2.4  and  ?.5l  are  intenraole 


[5],  and  admit  an  infinite  number  of  additional  differential 
conservation  laws.  For  )p>0  the  corresponding  equations  are  no 
longer  conservation  laws,  but  include  dissipative  terms. 


Solutions  for  =  0  are  found  by  the  method  of  1ST,  where 
a  certain  function  of  the  physical  amplitudes  is  used  as  a 
potential  in  an  associated  linear  wave  equation  (see  references 
given  above).  Two  formulations  of  the  initial  value  problem  are 
possible.  In  the  first  formulation  the  dependent  variables  are 
assumed  to  be  given  for  all  positions  at  some  initial  time,  and 
their  time  evolution  is  sought.  In  this  case  the  phonon 
amplitude  Q  appears  as  a  potential  in  the  associated  linear 
equation  [ 5 ] .  In  the  second  formulation  the  dependent  variables 
a~e  assumed  to  be  given  at  ail  times  for  some  initial  position, 
and  their  spatial  evolution  is  sought.  The  soliton  ootential  is 
now  a  function  of  the  optical  fields  [15].  While  this 
formulation  directly  solves  the  physical  initial  value  problem, 
where  the  optical  fields  are  given  at  the  entrance  to  the 
scattering  medium,  the  corresponding  conservation  laws  are  more 
difficult  to  handle,  since  they  involve  non  local  functions  of 
the  fields.  We  shall  therefore  use  the  conservation  laws 
obtained  from  the  first  formulation.  The  first  two  conservation 
laws  involving  the  pr.onor.  amplitude  Q  are  given  by  [5,24]: 

;  Q  q'it  -  C.5  S37  =  -  2  f  Q  Q*  ,  (2.6) 

Im(  Q*  Q ^  )x+  Im(  0. 5  Q* )  =  -  2  f  Im(  Q*  Q  ^  )  .  (2.?: 

The  first  equation  is  a  conservation  law  for  polarization 
energy,  balancing  the  time  change  of  polarization  energy  density 
in  the  medium  with  the  divergence  of  the  energy  current  carried 
by  the  optical  fields.  For  ^p=0  the  pump  is  depleted  in  space 
time  regions  of  increasing  medium  polarization,  and  depletion  is 
reversed  in  region  of  decreasing  polarization.  For  }f>0  energy 
is  lest  to  the  medium,  and  the  direction  of  energy  flow  for  the 
fields  depends  on  the  net  balance  between  the  rates  of 
polarization  change  and  energy  loss.  In  steady  state  only  the 
loss  ‘err.  is  kept,  and  the  pump  is  always  depleted. 

The  second  equation  gives  a  non  trivial  result  only  for  the 
case,  where  there  is  a  space  and  time  dependent  phase  difference 
between  fields  and  medium  excitation,  for  example  when  the 
optical  fields  are  not  exactly  in  resonance  with  the  Raman 
"ransition.  In  this  ‘ar,e  energy  and  momentum  nalar.ee  is 
maintained  through  the  generation,  of  an  op‘:ca,  pnor.on  wave. 
Equation  [2.1)  is  a  momentum  conservation  law  for  this  wave, 
which  balances  the  time  change  of  the  momentum  density  of  the 
phonon  wave  (first  term)  with  the  divergence  of  the  excess 
momentum  current  carried  by  the  fields  (second  term; .  The  latter 
current  is  obtained  from  12.2)  as 


In  the  presence  of  coherence  decay  momentum  is  dissipated  in  the 
medium,  which  is  described  by  the  right  hand  term. 


We  now  consider  the  one  soliton  solutions  for  the  case  r-  =o 

[16]: 


Q  =  exp (a B)  sech(A) 

Ap  =  if/,.  UR'  exp(:3 

As  =  V wR'  :  U 

S+  =  -  2  exp(iB) 

S3  =  (  2  §  sechpA)  - 

4  -  -  rif  ■ 

A:  -  “R'(  “r  *  u 


sech(A) 

tar.h !  A)  -  i  CJ  T  ] 
sech(A)  [  tar.h  (A) 

) 

3  =  O.T  - 

1  ,  =  O,  ;  co 


(2.9) 

(2.10 ) 

-  i  C*>T  ]  ,  (2.11! 


This  solution  describes  a  coherent  excitation  of  both  medium 
and  fields  travelling  at  a  speed  v  smaller  than  the  speed  c 
of  light.  Except  for  a  phase  factor  the  phonon  amplitude  is 
symmetric  about  the  temporal  center  (A=0).  Hence  gain  and  loss 
for  the  optical  fields  are  exactly  balanced  (see  2.6)  and 
photons  are  transferred  back  to  the  pump  pulse  in  the  trailing 
edge  of  the  soliton.  The  symmetric  shape  of  the  phonon  amplitude 
is  in  turn  caused  by  a  rapid  phase  shift  in  the  Stokes  field  at 
the  soliton  center.  The  excitation  can  be  observed 
experimentally  as  a  localized  pulse  of  pump  radiation  or  as  ar. 
(infinitely)  extended  pulse  of  Stokes  radiation  with  a  localized 
dip  in  intens i ty . The  temporal  width  of  the  excitation  is  equal  to 
IV  =  1  On  which  defines  an  intrinsic  time  scale.  The  solution 
will  be  a  valid  approximat  ion  to  (2.1)  and  (2.2)  if  jf"  <  <■  . 

The  frequency  of  the  pump  field  is  detuned  from  exact  Panar. 
resonance  by  &C0  =  .  I  f  Or  f  0  the  optical  fields  carry  a 


phase  wave  with 
soliton  center 

local  f 
is  equal 
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At  maximal  height  the  pump  pulse  reacr.es  a  fraction  §  of  total 
intensity.  Equations  (2.14)  summarize  the  observable  parameters 
of  the  soliton  and  their  relations: 


At  -  :  O  R  -  At o  =  to  z 
5  =  Mi  u p  =  :  : ;  1  *  2 at1  )  , 


1  'V 


,  jl  :  i  \ 

=  i/c  +  i/C toR  +  to t  ; 

The  three  components  of  the  optical  field  Bloch  vector  and 
the  local  frequency  are  shown  in  figure  1.  For  sufficiently 
large  negative  values  of  A  the  solution  closely  approximates 
the  physical  boundary  conditions  for  a  medium  extending  into  the 
positive  half  plane  %  >0  and  optical  fields  vanishing  or 
assuming  their  asymptotic  values  for  T<0  . 

In  the  presence  of  coherence  decay  solutions  exist  which 
show  similar  features  as  the  soliton  discussed  above.  In 
particular  localized  excitations  exist  which  show  reversal  of 
pump  depletion  ' 1 , 1 4 , 1 5 , 1 6 , 1 7 ’ .  If  the  typical  time  scale  of 
change  for  the  fields  is  larger  than  or  comparable  to  the 
coherence  time  (steady  state  or  transient  regime),  these 
solutions  may  appear  to  be  quite  different  from  the  soliton 
solutions  discussed  above  [25].  However  if  the  relative  phase 
between  the  optical  fields  varies  sufficiently  slowly  in  time, 
these  solutions  will  show  temporal  narrowing  and  eventually 
approximate  the  soliton  solutions  closely,  as  the  width  becomes 
much  smaller  than  the  coherence  time  (hypertransient  regime). 
Since  we  are  primarily  interested  in  the  dynamics  of  this 
approach  to  the  hypertransient  regime,  a  discussion  in  terms  of 
the  exact  soliton  solutions  for  the  solvable  case  is  appropriate 
and  useful . 

As  a  result  of  coherence  decay  temporal  width,  position  and 
maximal  amplitude  will  be  general  functions  of  propagation 
distance,  and  the  pulse  shapes  will  no  longer  be  symmetric.  It 
is  found  that  to  first  order  in  the  changes  in  width, 

position  ar.d  amplitude  can  be  calculated  without  taking  into 
account  the  changes  in  pulse  shape  [21].  We  shall  calculate 
these  parameters  by  integrating  equations  2 . 6  and  2.9  over  time 
This  gives  ordinary  differential  equations  in  *  for  the 
corresponding  integrals.  The  integrals  are  calculated  by  using 
solutions  2.9  and  2.11,  however  with  parameters  00%  and  COx 
that  are  functions  of  ^  •  The  resulting  equations  are: 

(d-d*  )  Ik  =  -  2  r  Jk  ,  k  =  1.2  ;  (2.15 
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The  solutions  are  most  easily  expressed  in  terms  of  the 
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The  coefficients  involving  the  propagation  distance  can 
be  expressed  in  terms  of  the  steady  state  gain  G  or  the 
transient  gain  G'  .  The  transient  gain  G'  is  a  measure  for  the 
spatial  change  of  the  soliton  fields  and  is  appropriate  for  a 
discussion  of  the  internal  dynamics  of  soliton  propagation.  It 
is  defined  from  equation  (2.12)  as  G'  =  2  .  On  the  other 

hand  for  soliton  experiments  the  steady  state  gain  G  is  more 
relevant,  since  its  value  is  limited  by  processes  like  higher 
order  Stokes  generation  to  (approximately)  G  <  25  .  This 
parameter  determines  the  growth  of  Stokes  intensity  for 
undepleted  pump  in  the  steady  state  regime: 

:s  =  Is(0)  exp(G)  ,  G  =  2  /  f  .  (2.19) 

In  terms  of  G  the  coefficients  in  (2.17)  and  (2.13)  are 
given  by: 

’  2  f2  G  ’  =  2  !  if  hi'2  =  ■  !2-20: 


For  exact  Raman  resonance 
width  from  (2.18) : 
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In  the  limiting  case  of  small  gain  and  width  the  width  is  seen  to 
decrease  linearly  with  gain.  The  corresponding  coefficient  is 
proportional  to  the  square  of  the  initial  width.  In  the  opposite 
limit  of  large  gain  the  width  becomes  independent  of  its  initial 


value,  and  decreases  with  the  inverse  square  root  of  the  gain. 
This  limit  is  not  actually  reached  for  realistic  gain  values  in 
the  hypertransient  regime.  Numerical  results  show  however  that 
the  soliton  width  does  indeed  become  almost  independent  of  the 
initial  width  at  large  gain. 

For  non  zero  detuning  the  soliton  amplitude  g  decreases 
exponentially,  while  the  parameter  Ja^  ,  whose  inverse  is  the 
maximal  local  frequency  of  the  corresponding  phase  wave,  remains 
constant.  If  the  soliton  width  is  of  the  order  of  the  coherence 
time,  and  the  detuning  is  much  smaller  than  the  line  width,  we 
have  : 

jtt.  R  £  AU  (AX)2  ,  L  *  2G(4(0f(AT)2)2;  (2.22) 

showing  that  the  attenuation  coefficient  L  for  the  soliton 
amplitude  is  proportional  to  the  square  of  the  detuning.  This 
explains  the  observed  rapid  decay  of  Raman  solitons  for  off 
resonant  SRS  [14].  The  soliton  width  reaches  its  minimum  for 
^  =  0.5  and  increases  after  that  point. 


III.  SOLITON  PROPAGATION  IN  I NHOMOGENEOUSLY  BROADENED  MEDIA 


In  inhomogeneous iy  broadened  media  the  matrix  element  Q  , 
which  determines  the  medium  polarization  at  the  optical 
frequencies  is  obtained  as  a  statistical  average  over  molecular 
subpopulations,  each  with  different  frequency  shift  A  from 
exact  Raman  resonance  in  the  laboratory  frame.  The  statistical 
distribution  is  characterized  by  the  line  shape  function  g(A  ). 
If  we  denote  the  ma"rix  element  for  the  population  with  frequency 
shift  A  by  Q*  equation  2.4  for  the  matrix  element  is 
replaced  by: 


1 A  g  [  6 


(3.1) 

(3.2) 


We  integrate  (3.1;,  substitute  -  he  result  into  (3.2)  and  obtain: 
QA  (  zr  )  =  |  dr'  exr- '  i  A  I  r  -  t  '  )  1  Q .  5  S  (  t  ')  ,  (3.3) 
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g(  t  )  =  j  d  t  exp [  i  &  t  ]  g( A  ) 


line  shape  function  g 


is  the  Fourier  transform  of  the  spectral 
(temporal  correlation  function). 


We  note  that  the  case  of  homogeneous  broadening  is  included  in 
(3.4)  with  g{  x  )  =  exp(-J-T  ).  This  shows  that  in  the 
absence  of  medium  saturation  homogeneous  broadening  is  equivalent 
to  inhomogeneous  broadening  with  Lorentzian  line  shape.  The 


conservation  laws 
of  Q  ,  which  is 


(2.6)  and  (2.7)  require  the  time  derivative 
obtained  from  (3.4)  as: 


j  d  r'  g  (  t  -  x' 


3.5) 


where  g  is  the  derivative  of  g  .  The  last  term  in  (3.5)  is 
the  correction  from  coherence  decay  to  the  exactly  solvable  case. 
We  shall  discuss  here  only  the  case  of  exact  resonance  at  line 
center,  and  of  a  symmetric  line  shape.  In  this  case  Q  is 
real,  and  we  obtain  from  the  first  conservation  law  in  analogy 
with  (2.6): 


Q  Q). 


(3.6) 


The  correction  term  is  calculated  to  leading  order  in  the  line 
width  by  expanding  the  correlation  function  g 


=  g ( o )  Qv  '  +  g ( o ; 


Q,0,(r') 
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Here  q(°)  is  the  matrix  element  for  vanishing  line  width.  For 
homogeneous  broadening  the  first  term  gives  a  non  vanishing 
result:  g(C)=  -  J-  ,  due  to  the  weak  decrease  of  the  Lorentzian 
spectrum  at  large  frequencies.  For  inhomogeneous  broadening  on 
the  other  hand  with  a  Gaussian  line  shape  or  other  shape  with 
sufficiently  fast  decrease  at  large  frequencies  the  first  term 
vanishes,  and  the  leading  contribution  is  of  second  order  in  the 
line  width.  In  this  case  we  obtain  from  (3.6)  the  differential 
equation : 
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For  the  case  of  a  rectangular  spectrum  with  half  width  the 
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The  relation  between  propagation  distance  and  steady  state  gain 
is  in  this  case : 


g(0)^ 
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giving  the  final  result  for  the  soliton  width  as  a  function  of 
gain : 

AT  =  At  (0)  (  1  -  0.5  -rr  (  ^AT(0))3G)(_1/3)  (3.12) 

AI(O)  (  1  -  (TT/6)  (  f  AT(0))3  G  )  . 

The  rate  of  decrease  is  now  proportional  to  the  third  power  of 
the  initial  width,  and  the  small  gain  limit  considered  in  the 
second  line  of  (3.12)  above  is  the  only  realistic  case. 

We  have  not  studied  the  off  resonant  case  so  far.  Numerical 
results  indicate,  that  the  behavior  can  be  quite  complex  in  this 
situation.  For  sufficiently  large  detuning  width  and  amplitude 
can  decrease  and  increase  alternatively,  and  strong  pulse 
deformations  may  develop.  We  shall  discuss  these  results  in  the 
following  section  IV. 


IV.  NUMERICAL  RESULTS 

The  transient  SRS  equations  (2.4)  and  (2.5)  were  solved 
numerically  with  the  one  soliton  form  (2.11)  for  the  optical 
fields  as  functions  of  time  as  initial  condition  at  zero 
propagation  distance.  The  soliton  amplitude  and  width  were 
determined  for  the  numerical  solutions  from  the  maximal  value  of 
S  and  from  the  integral  I  ( equation ( 2 . 16 ) .  For  inhomogeneous 
broadening  a  rectangular  line  shape  was  chosen  for  numerical 
convenience,  and  the  line  was  sampled  at  a  finite  number  of 
equally  spaced  points. 

Figure  2  compares  numerical  (circles)  and  analytical  (solid 
lines)  results  in  the  hyper  transient  regime  and  the  resonant  case 
for  homogeneous  (lower  curves)  and  inhomogeneous  broadening 
(upper  curves).  The  agreement  is  excellent  in  the  latter  case, 
and  better  than  10%  for  the  homogeneous  case.  In  this  and  all 
subsequent  cases  the  time  scale  is  fixed  by  choosing  the  initial 
soliton  width  as  the  unit  of  time.  The  line  width  is  0.:  for  the 
homogeneous  and  0.1571  for  the  inhomogeneous  case,  which  results 
in  the  same  gain  per  propagation  distance  and  the  same  coherence 
time  for  both  cases.  The  soliton  width  is  hence  only  10%  of  the 
coherence  time.  At  a  gain  of  20  the  homogeneous  soliton  has 
narrowed  by  about  15%,  while  the  inhomogeneous  soliton  has 
narrowed  only  by  about  3%. 

Figures  3  and  4  compare  numerical  (circles)  and  analytical 
(solid  lines)  results  for  soliton  width  and  relative  amplitude  in 
the  off  resonant  case  and  for  homogeneous  broadening.  Initial 
soliton  width  and  line  width  are  as  before  and  correspond  to  the 
hypertrans lent  regime.  The  values  chosen  for  detuning  are  0.1, 
0.5  and  1.1.  These  values  correspond  to  the  lower,  middle  and 
upper  set  of  curves  in  figure  3,  and  to  the  upper,  middle  and 


lower  set  in  figure  4.  As  predicted  the  soliton  width  (figure  3} 
is  larger  for  larger  detuning.  For  the  first  two  cases  is 
actually  larger  than  predicted,  with  a  predictive  accuracy  of 
about  10%  .  Somewhat  larger  discrepancies  occur  for  a  detuning  of 
1.1  .  Here  the  soliton  width  is  predicted  to  increase,  while 
the  actual  width  almost  remains  constant.  A  similar  behavior  is 
found  for  the  soliton  amplitude  (figure  4).  Larger  discrepancies 
occur  again  for  a  detuning  of  1.1  ,  where  the  actual  amplitude  is 
larger  than  predicted.  These  discrepancies  are  caused  by  the 
behavior  of  the  phonon  wave.  In  particular  the  local  frequency 
(see  equation  (2.13))  begins  to  show  strong  deviations  from  the 
one  soliton  form  (2.13)  at  larger  detunings,  including  an  ai.mos~ 
linear  frequency  chirp  with  increased  detuning  at  the  trailing 
edge  of  the  soliton. 

The  effects  of  detuning  are  much  .more  pronounced  for 
inhomogeneous  broadening,  especially  if  the  detuning  is 
comparable  to  or  exceeds  the  line  width.  This  is  plausible  in 
view  of  the  fact  that  the  weight  of  higher  frequency  componerhs 
is  much  smaller  in  the  inhomogeneous  case  than  in  the  homogeneous 
case.  A  thorough  physical  analysis  and  analytical  results  remain 
to  be  given.  Figure  5  shows  numerical  results  for  sol  item  width 
(solid  lines)  and  amplitude  (dotted  lines)  for  a  rectangular  line 
with  a  half  width  of  0.1571  and  an  initial  soliton  width,  of  1.0 
( hypertransient  regime).  The  values  for  detuning  are  C.Zi,  C.l 
and  0.2.  At  given  propagation  distance  the  soliton  width 
increases  and  the  amplitude  decreases  with  increasing  detuning. 
For  small  detuning  (0.05)  the  width  decreases  as  a  function  of 
distance,  while  it  increases  initially  for  large  detuning  (0.2) 
and  has  a  maximum  at  a  gain  of  about  12.  At  about  the  same  gain 
the  soliton  amplitude  reaches  a  minimum.  This  behavior  is  in 
strong  contrast  to  the  homogeneous  case. 

Soliton  narrowing  is  much  stronger  in  the  transient  regime, 
where  the  temporal  soliton  width  is  comparable  to  the  coherence 
time.  In  the  analytical  treatment  above  we  found  that  the 
narrowing  rate  per  gain  unit  is  proportional  to  the  second  power 
of  the  line  width  for  homogeneous  anc  to  t r.e  third  power  for 
inhomogeneous  broadening.  Figures  6  and  7  show  numerical  results 
for  width  and  amplitude  for  homogeneous  anc  inhomogeneous 
broadening.  The  initial  so.itor.  width  is  l.C  and  the  line  width 
is  1.0  for  homogeneous  and  1.571  for  inhomogeneous  broadening. 

The  values  for  detuning  are  3. 05,  0.1  and  0.2. 


Tor  homogeneous  broadening  (figure  6)  a  detuning  of  0.05  or 
5%  of  the  line  width  must  be  considered  a  limit  for  practical 
app  1 1  ca  t  .or.s  to  pulse  narrowing.  In.  this  case  the  soliton  decays 
to  about  90%  at  a  gain  of  20  and  narrows  to  a  temporal  width  of 
20%  of  its  in.it. a*  va.ue.  At  10%  detuning  the  90%  amplitude 
level  is  reached  already  at  a  gain  of  7,  and  the  soliton.  wid*h 
is  40%  of  tr.e  ir.itia.  value.  If  is  interesting  to  note  that  th<= 
initial  rate  of  soliton  decay  is  much  less  than  in.  the 


to  note 
in.  the 


ha  t 


hypertransient  regime,  where  exponential  decay  is  found. 

Inhomogeneous  broadening  figure  7  leads  -.o  better 
stability  for  solitor.s  at  sma.l  detuning  than  homogeneous 
broadening.  At  a  detuning  of  C . C5  or  1%  hie  soli  tor  amplitude 
reaches  a  minimum  f  9T%  a-  a  gain  uf  12.  »r.d  .c  reases  after 
that  point.  The  nar-iwm.g  is  a.  mis-  tne  same  as  m  me 
homogeneous  case,  w  1  - a  larger  rate  of  narrowing  at  sma.l  ga  :  n  . 
At  a  detuning  of  C.l  or  6%  »he  so .  ;  tor.  amplitude  stabilizes  a" 

70%  For  larger  det  ur.  i  r.g  7.7  r  1  ’’%  an  .  n.  termed  late  man. mum 
in.  amplitude  occurs  at  a  ga.t  jf  .  v  and  n.e  wi  inn  begins  to 
increase  rapidly.  A*  *  -  .  s  :•■."*  .  -ng  range  osc  i  ,  .  a  t  i  or.s  begin  to 

develop  in  the  ‘rai.i  r. g  euue  I  *  ne  so  .  i  -  m.  . 

Tempora.  pu.se  sr.apes  for  the  pump  beam  are  shown  i  r.  fig 
3  and  9.  Figure  o  snows  n  rrogereous  or  Jdien mg  for  a  defun m 
9.0  !  0  .  a  :  ,  C  .  Oh  '  6  .  o  ;  .  7.1  3  .  c  ;  ana  C  .  2  3  .  d  ;  .  Figure  9  sr.v 

inhomogeneous  nr  .dien.t;  fir  a  detuning  if  C . 0  9.  a1,  C.15  9 

0.1  (  9  .  c  )  and  T  .  2  9  .  u  .  Strong  ,sc  ;  .  lat  ions  are  seen,  to  j.  - 

in  the  framing  edge  fir  inhomogeneous  oroader.  i  r.g .  These  fade 
away  at  .a rger  gain  f  r  zeru  and  sma  1  i  detuning  ! 9 . a  ana  9 . b i 
but  persist  for  .arger  ue-uning  9 . r  arc  9 . d )  . 

V.  S'JMMAPV 

Pamar  soli  tons  correspond  exper  i  men  t  a  .  i  y  *o  short  localized 
pu.ses  of  pump  rename:.  in  an  envelope  of  Stokes  radiation. 

They  are  sustained  by  Da  lamed  processes  c  f  absorbt  ion  and 
emission.  In  the  preserve  of  coherence  de^ay  pulse  narrowing  and 
decay  occur.  Dec  ay  *  urs  only  if  the  -wo  apt  ir.nl  fields  are  net 

ex.ac-. y  =r.  p4axar>  rescr.ar.  :e  '"detuning";  . 

We  nav  ,*  .!.>■-  :  *  >-se  •  r.e-  me -a  n  *  .  n..  by*  leal  ly  and 

numerically  .  r.  the  v/perf  rats.-rv  regime,  a n.c  tuner  i-:a  1  ly  .:.  tne 
‘  vans  lent  regime  .  The  ases  of  '00*;-  me  gene;  .us  and  1  nhemogeneous 
oroader.  mg  w-re  s*  ud.ei 

'ur  in.a  v;t  vai  *  ^  vcp-*  wn  need  ~  ~  >r.t  :  a  1 

.  .r.at'Vi*  1  n  .-.wo  f  or  p ro'  er.e  -gy  a~d  .r.o.m.er  *  um  .  Ir.ese 
conservation  .  aws  g.ve  mat  r.*T.a*  .  ».  expression  *  •.  t  he  ba  .am  e  of 
pr.ysical  processes  an  i  *h*»  effe ■■  t ,f  nroacer.  ,:.j  :.  m.is  ba. ar.ee. 

For  the  off  resonant  ase  f  -«»<p;en-  y  detun.  nj  *  he  occur  -er.ee  of 
a:,  jptical  phonon  w.v-  .  s  m  maintain  ;y  and  r..,;r.er.fur. 

conservation,  ana  *0  have  ;  >  vis  .  de  rao  1  e  mi.  ie:  e  sa.  item. 

3  d  C  I  -  1  t  V  . 

In  t~e  t  rar.s  1  erv  'eg. me  we  f.rd  a  ie*  1  j  f  1  %  ;  n  1  •  he 

Parnan  ..n.e  w.ctrv  *  0  r-ys^nt  a  praC.cal  1  .m.v  for  tne 
acr.ievement  of  m  apt  r>-  1  a  n  .  e  na-r  w  .  ng  and  g-  .  1  w.  1  i  ♦  on 

stability.  Solitor.s  are  mere  stab .  e  for  inhomogeneous 
broadening,  however  r>s<  .  imiuv,  occur  in  -r.e  ■  ;  -rg  edge  of 

me  pump  pulse  and  persist  for  .anger  detunings. 
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the  off  resonant  case  are  shown  for  inhomogeneous  broader.! 
the  transient  regime.  The  Raman  line  width  is  1.571  and  t 
detuning  is  0.05,  0.1  and  0.2.  In  this  sequence  of  values 
lower,  middle  and  upper  solid  lines  are  the  width,  while 
upper,  middle  and  lower  chain  dotted  lines  are  the  soliton 
amplitude.  Frequencies  are  in  units  radians  per  time  unit, 
the  time  unit  is  the  initial  soliton  width. 

Figures  8a  and  8b  The  temporal  pump  pulse  for  homogeneous 
broadening  in  the  transient  regime  is  shown  for  a  detuning 

;8a)  and  0.05  (3b).  Pulses  a-e  shown  for  a  gain  increment 

beginning  with  the  initial  pulse  in  the  upper  right  part. 

Figures  be  ana  3d  1  r.e  tempo ra*  pump  pulse  for  homogeneous 
broadening  in  the  transient  regime  is  shown  for  a  detuning 
(8c;  and  C.2  (3d).  Pulses  are  shown  for  a  gain  increment 

beginning  wi*n  the  initial  pulse  in  the  upper  right  part. 

Figures  9a  ana  9o  Tr.e  temporal  pump  pulse  for  homogeneous 
broadening  in  the  transient  regime  is  shown  for  a  detuning 

9a;  and  C.C5  ;9b) .  Pulses  are  shown  for  a  gain  increment 

beginning  with  the  i  r.  *  *  i  a  1  pulse  in  the  upper  right  part. 

Figures  9c  and  9d  The  temporal  pump  pulse  for  homogeneous 
broadening  ;  r.  the  ‘•ransierd  regime  is  shown  for  a  detuning 
?c  and  0.2  :3d'.  Pulses  are  shown  for  a  gain  increment 

beginning  with  the  initial  pulse  in  the  upper  right  part. 
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APPENDIX  C 


Effect  of  Medium  Saturation  and  Coherence  Relaxation 
on  the  Propagation  of  Raman  Soli  tons 


Christian  J.  Tourenne  and  Kai  J.  bruhl 
Department  of  Physics 
slaharishi  International  University 
Fair tie id,  Iowa  5255b 


ABSTRACT 

The  effects  of  medium  saturation  and  homogeneous  broadening  rn 
the  propagation  of  soli  tons  in  stimulated  Raman  scattering  are 
discus seil .  The  general  problem  including  Stark  shift  and  medium 
saturation,  but  without  broadening  has  U-en  solved  by  Kaup  and  Stru¬ 
del  using  an  extension  of  the  AKNS  methcxi.  For  the  limiting  ease  or 
negligible  Stark  shift,  we  obtain  first  enter  equations  for  the 
sol  it  on  width  in  the  presence  of  homogeneous  nroadening.  We  use  t  uese 
asymptot  ie  solutions  to  discuss  the  effect  of  metium  saturation  on 
pulse  narrowing  and  eompare  them  to  a  numerical  analysis. 


In  stimulated  Raman  scattering  (SRS )  energy  is  exchanged  between  two  laser 
beams  at  different  frequencies  bv  nonlinear  interaction  in  a  suitable  medium. 
Energy  is  normally  transferred  from  the  high  frequency  Lieam  (pump)  to  the  low 
frequency  beam  (stokes).  In  both  numerical  [1]  and  laboratory  [  ‘j. ,  3 1 
experiments,  so  1 i ton  1 i ke  excitations  have  been  created  bv  a  sudden  phase  change 
of  1H0  degrees  in  t.he  injected  stokes  beam .  In  the  presence  of  coiiisionai 
coherence  decay  in  the  medium  and  for  exact  resonance,  the  pulse  shows  t  no 
features  of  temporal  narrowing  and  acceleration  in  the  laboratory,  frame 
[1.3,41. 

The  effect  of  coherence  decay  on  the  propagat l on  of  Raman  so i i tons  has 
already  been  stud i ed  txith  analytically  and  numerically  using  asymptotic- 
perturbation  theory  [i,3,lj.  In  these  studies  however,  certain  effects  like 
medium  saturation  or  dynamic  Stark  shift  were  neglected. 

In  this  paper,  we  discuss  how  medium  saturation  affects  temporal  puls* 
narrowing,  using  a  perturbative  method  which  is  based  on  the  energy 
conservat ion  law  for  medium  and  optical  fields  [5]. 

The  general  problem  of  soliton  propagation  in  stimulated  Raman  scattering 
(SRS),  including  Stark  shift  and  meditim  saturation,  has  been  solved  by  kaup  [61 
and  Steudel  [7]  using  an  extension  of  the  \KNS  method  [8|.  In  the  presence  of 
coiiisionai  coherence  decay  and  for  the  limiting  case  of  negligible  Stark 
shift,  kaup’s  and  Steudel ’s  equations  [6,7]  for  transient  SRS  become: 
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t  H'.ils,  r  ;s  profx.rt  lonal  to  t  lie  population  d  1 Merer*  •**  between  t  in*  «  wo  i 
f  trie  Raman  transition,  sj  :s  proport  l  ona  1  to  t  tie  iiffernce  of  ptiot  >n  -urr*  'tit 
densities  net  ween  the  pump  and  t  he  stoius  tx-ams  uiu  s  is  nrop-ort  lonai  to  •  he 
prodix  t  of  the  stokes  elec-trie  fedd  ;uid  the  complex  ,'i  m.lliCit  e  of  ’Me  pump 
"let  ric  field.  V  tnd  \  ire  'ime  .uni  space  like  c<x>r  iinat  es  m  a  re-  mied  »  .me 
t  rnme  and  t  lie  subscript:-  indrate  ear'  lal  di  t  ferent  iat  ion  with  res|»vt  to  nil's 
'  oor-r )  i  nat  es  .  i  ni.-.euent  unit-,  nave  tx-en  isosen  so  'hat  all  oupl:ng 
•  •<  itist  ant  s  at-"  equal  ;  me  tn.  *  erm  -  K  t  desi  r  i  i*>s  ••  d  i  i  s  i  >na  ,  .  •<  >h<  ren< -e 

(e,  wi!n  eoherenie  tune  ■  y  ,  chef.  J  is  the  imlui.u  :  .iman  ;ir,e  wid’h.  m 
have  not  onsidored  jx.pui  C  .  ,i;  r< jaxit  ion  i.ete,  uni  ....  is  Minn • 

'  .Hie  sea  1  e  inuefi  :  .tiger  trial:  till  oherence  *  :  me  . 


I  I 't 


In  the  hypertransient  regime  (<r+)  >>Jr+),  the  phvsical  quantities  vary 
on  a  time  scale  much  shorter  than  the  dephasing  time  1/y  ,  In  the  iimit  y  =  u, 
equations  (1.1)  to  (1.4)  are  integrable  and  exact  solutions  can  be  found  using 
an  extension  of  the  inverse  spectral  transform  [6,7).  For  the  special  case  of 
purely  imaginary  eigenvalues,  the  one-soliton  solutions  to  equations  (1.1)  to 
(1.4)  are : 


2  H  V  2 

-(f}j  -1/4)  cosh<2j*)/(  (Q1--l/4)cosh  (2 

-1  +  (  l/2)/i  2-l/4)cosh2(2|i)  +  l/4) 

(  2n,  l  n  ,  ^—1/4  ?^sinh(  2 it)  )  /  (  (  rt ,  ^-1  /4  )c< 


jjJ  +  l/'4> 


2^1^  -l/4f»sinh(2i».)  >/(  (^  -1/4  )cosh"(  2u.)  +  l/4  ) 


-1  +  2y^/(  i  ^,  “-1/4  )cosh^(  2^*,)  +  l/4  ) 


(2.1) 


(2.2) 


(2.2) 


(2.4) 


where 


is  the  imaginary  part  of  the  eigenvalue  and  -  X/4IJ .  . 


In  this  case,  the  optical  fields  are  exactly  in  resonance  with  the  Raman 
transition.  In  the  foLiowing  we  discuss  the  physical  processes  underlying 
soliton  propagation  for  different  values  of  ^  .  in  the  leading  edge  of  the 

pump  pulse  (  |J->0),  s  +  is  in  phase  with  the  induced  polarization  r+  in  the 
medium,  leading  to  a  decrease  in  s^  (stokes  scattering).  In  the  trailing  edge 
(J*-<0),  however,  s+  has  changed  sign  and  s.^  increases  now  with  propagation 
distance.  This  means  that  the  stokes  beam  experiences  loss  with  corresponding 
gain  for  the  pump  beam  (anti  stokes  scattering).  By  the  same  mechanism  t  tie 
medium  is  excited  in  the  leading  edge  ana  deexciten  in  the  trailing  edge.  Both 
processes  (loss  and  gain,  excitation  and  deexcitation)  are  perfectly  balanced, 
as  is  evident  from  the  symmetric  soli  ton  shajx*. 


Two  limiting  cases  are 


i nt  erest  : 


a)  when  ,  becomes  large  (^.  -+«o) ,  t  he  almve  anaivt.icai  solutions  ire 
equivalent  to  t  he  one  soliton  solutions  found  Bv  Chu  and  Scot  t  j  a  )  for  t  tie  ease 
in  uhieh  medium  saturation  is  neglectful  ( r  — ♦  -1);  m  this  case,  t  he  temjjural 
width  of  the  so  i  1 1  on  is  mueh  smai  ier  than  the  Rabi  t  mie  ,  during  wh  e  n 
|x>pulation  inversion  ixrurs,  and  the  gain-ioss  reversal  net  ween  the  -.tones  anti 
the  pump  teams  occurs  very  rapid iy; 

b)  when  rj  approaches  1/2,  t  ne  rise  i  une  of  s  i*  forties  comjxirahl  >■  to  tin' 
kati  i  time  leading  to  a  rxipulat  ion  inversion  for  the  medium.  The  of  f-di  agonal 
matrix  element  r  is  cor  resound l ng 1 \  small,  leading  t  ■  a  more  gradual  change  it 
the  t  leld  variables.  As  a  eonsequerice .  a  wide  temporal  plateau  i  t  ortur'd  . 


:n  the  transient  regime  Mr  no  analvtieal  solution  of  equations 

!.!)  ti.  ti.l)  exists  in  <  •  1  used  form  i  e\eept  m  t  )ie  linear  regime  |iOi,  where 
iteplet  i  on  of  the  pump  lie  am  is.  neg  1  i  g  1 1 1 1  e  and  the  pump  amplitude  is  assume!  to 
lie  mdeix-ndent  of  x).  However,  for  a  finite  damping  rate  (  >0),  the  pnvsieal 

mef'hanism  resjixinsi  ble  toi  soliton  proiiagat  i  on  m  the  hypert  rtins  l  ent  regime  is 


expected  to  be  valid  also  for  transient  SRS  (  any  phase  change  for  s  will  lead 
to  gain-loss  reversal  and  solitonlike  excitations). 

If  the  coherence  decay  time  T  =  1/y  is  long  compared  to  the  pujse  width, 
approximate  solutions  of  equations  (1.1)  to  (1.4)  can  be  obtained  through 
perturbation  theoretic  methods  fill.  Under  experimental  conditions  [lj,  the 
optical  fields  are  given  as  functions  of  time  and  one  wants  to  find  their 
evolution  in  space.  Kaup’s  general  perturbation  theory  [ 1 1 J  cannot  therefore  be 
directly  applied,  since  the  soli  ton  is  defined  in  terms  of  the  optical  fields, 
while  the  perturbation  occurs  in  the  equations  for  the  medium  variables.  A 
generalization  to  this  case  seems  possible,  and  has  been  given  for  the  limiting 
case  a)  above  by  kaup  [12].  We  shall  here  use  a  more  simple  and  direct  approach 
which  uses  the  following  energy  conservation  law  derived  from  equations  (1.2) 
and  (1.4): 


(S3 


-  (r3>T 


(3) 


The  basic  assumption  of  our  method  is  that 
solutions  of  equations  (1.1)  to  (1.4)  have 


for  smal  1  , 

the  form: 


the  analytical 


r  =  r° ( n , ( X ) )  +  r 1  (4.1) 

r3  =  r3  (l7l 1  *  *  1  +  r3  (4.2! 

s+  =  s«in  (X)>  +  SJ  (4.2) 

s^  =  s*( ^ ( X ) )  +  s *  (  4.4) 


The  first  (zeroth  order)  terms  are  the  unperturbed  solutions  (2.1)  t.o  (2.  s)  in 
which  ^  is  assumed  to  be  a  function  of  X,  and  the  second  terms  in  the  above 
equations  are  first  order  corrections  to  the  one  so  I i t  on  solutions,  which 
deserilie  change's  in  shape. 


By  integrating  t  he  conservation  iaw  (3)  1**  tween  !'=-  *0  and  T=+«0  ,  w 
>bt  ai  n : 


ld/riX) 


-00 


r.  ( +«0  )  -  1 


i 


I'he  t  .me  integral  of  s  represents  trie  temporal  w.  It  h  w  •('  the  optc-i!  pulse: 


rJ-  oo 

=LS:> 


ciT 


In  onler  to  calculate-  the-  asymptotic  value  i  f  r.'+cOi,  we  w-,e 
( 1.2)  to  i  >btain : 


and 


I  r  *" )...  +  ■  r  “  i  ..  -  -  2  if  t 
+  i  a  1  + 

1  rit  egrat  i  ng  both  sides  from  F'-«0  to  :  =  +  00  .  end  r-  mo  ■mi»-r  i 1  -.-1  t  hat  *  he  >  t  *  — 

diagonal  mat  r  ix  eiement  should  vanish  at  T-  +  <®  "i  -  eO  .  we  g.-t 
equat  ion : 


'he  t  o  .  I  ow  l  ng 


r^(+00)  -  1  =  -  it  I  r+^  dT  (t  ) 

oo 

since  r^(-00  )=-l.  Equation  (6)  shows  that  the  asymptot  ic  valu**  .it'  r.^  at  T  =  +oo 
is  not  equal  to  -1.  Due  to  coherence  decay,  the  perfect  balance  between  medium 
excitation  and  deexcitation  is  upset:  at  the  trailing  edge  of  the  pulse,  the 
medium  polarization  is  smaller  than  at  the  leading  (dge.  As  a  result  ,  .a 
fraction  of  atoms  gets  trapped  m  the  upper  level  of  the  Raman  transit  ion. 

To  first  order  in  if ,  the  population  difference  r  ,+ooi  is  given  bv : 


r.j(  -too  )  -  1  + 


'S—  can 


Now  combining  equations  <  :>  >  and  (7),  we  obtain  our  final  result  fur  the  sjiat  :  a 
dependence  of  the  soli  ton  width: 


dw/oX  =  -  "S 


r  dT 


We  require  that  the  pulse  width  w  is  entirely  det ermined  bv  the  one  s  i i t  .n 
part  s  of  the  solution:  — , 

w  -  J  «;  .it 

w/_  oo 

To  first  order  the  right  hand  side  .  t'  equation  (Hi  is  det  ei-mi  ned  by  the 

one  soli  ton  contribution  onlv;  hence : 

,+  ° O 

dw/d\  -  2T J  r®“  dT  hi 

oo 

in  which  tr^  is  now  •»  function  of  X.  The  two  integrals  are  ,  -as  l  i  caici.»at»*d  tr. 
using  the  solutions  [2.1)  and  .11.1)  and  t  he  f.  1  lots  mg  restii  t  s  are  btamed : 


dw/dX  -  (d/d X)/  (s*+l)  dT  =  d/dX  |JiJi(i2r^ 


+  1  i  /  :  ;.o,  - ;  i  i  j 


r*Z  dT 

+ 


l  l/kJfJj'  )  |  l  +  ( 


?i+i  r 


i'he  temfxiral  width  of  the  sol  i  ton  is  related  i  <>  the  itaagiuarv  part  rj  i  ’  n< 

ej  gen\alue  by:  u=2Ln|  (2*^  +  1  )/(2Qj-l  1  |  .  Kquat  n»n  t  b  i  t  <»ge*  per  u  :  t  n  1  <*■  uio 
ill),  lead  to  the  following  differential  equation  for  write.  w: 


dw/dX  =  4  Jf<  exp(  w ) awexpi w/2 ) -  1 ) ( 2e\pt w/U > -cxpl w i - ! 


I  ’  l  e\p.  W : *_o\p 


This  equation  has  two  interesting  limiting  cases: 

*  when  H  is  large  (negligible  medium  sa'urat;  n ) .  : a. 

inversely  projjort.  i  ona  1  to  'he  square  root  of  X: 


w(  X  )  2/  (  X  a  1/w  [())-) 


in  accordance  with  a  previous  perturbation  analysis  of  the  unsaturated  case 

[91; 


*  when  the  medium  has  a  large  medium  saturation  (  — •  1/2),  the  soliton 

width  decreases  proportionally  to  X: 

w  ( X )  as  -  4f\  +  w  (  0  )  ( »|  j  — *  1  /  2  )  (14) 

In  order  to  check  the  accuracy  of  our  perturbation  method,  we  have  solved 
equations  (1.1)  to  (1.4)  numerically;  the  Raman  medium  was  assumed  to  be 
located  in  the  half  space  \>0  and  not  excited  initially;  for  optical  pulses  f 
finite  duration,  the  boundary  conditions  were: 


for 

T>0 : 

r  <X,T=0> 

=  i) 

r3(X,T-0) 

-  0 

s  ( X-0 , T ) 

+ 

=  s  (X=0 

+ 

s  (X=(),T) 

=  s  (X=() 

for 

T<0 : 

s  ( X -0 « T ) 

:  0 

s,  (  \— 0  *  T ) 

=  1 1 

['he  initial  fields  were  taken  to  ne  equal  to  their  >ne  soi  it  on  form  at  t  he 
ent  ranee  of  the  Raman  cel!  i  \=0  i  ano  their  evolution  was  then  <  -a  leu :  at.ed ,  is 
the  pulses  propagated  into  the  cell  t  X/0 ) ,  by  soi\  mg  (1.1)  :  o  (1.1) 
numer  i '  a  1  1  >  ,  for  different,  initial  values  of  and  for  different  sma  i  i  va  ues 

of  *  ne  homogeneous  t  To.aden  i  nil  (laramet.er  ¥  .  The  so  li  •  ;  in  widtti  given  ov  t  he 
numer'ioai  solutions  for  s  tor  different  values  T  the  gam  iengtn  \  was  t  hcti 
onumred  to  the  arialvii  aT  soiut  ion  for  the  ’  -  mj  x  r.i  i  widtti  given  .a  '-qua-  n 


!  or 
f  >1 . t  t 

.  at  t  IT 


fig.  .  mo  2,  we  shew  how  the  analvtii-ai  an<t  numerical  r*»st.i  t  •<  imuir*- 
[_  ,  TS  =u.l  and  rj  ro.  rrspivt  i  v  e  1  \  .  lie  w  idt  h  was  • ‘h<  >ser.  '  > 

sol  i  tot)  we.  1  within  t he  area  of  mtegrat  ior.,  v  cod  i ng  -r r  ■  t  -  t  r  ■  -rr.  •  ti< • 
r'oi  i  -di  jr  e  exjjonent  i  i !  I  v  sma  i  i  .  !'h«  •  soi  id  <  nirv  •  ano  the  j«».rit  .  mr  are 


the  numerical  and  ana  1  vt  i  <  -a  i  w  idt  hs  respe<  ’i.elv.  as  i  •  it  n  r  i  <  r  f  v  .  'or  *  n< 

'  •  -it. ;  k  >ra  1  kaman  so  1  ;  t ,  ms  .  At  h  i  ghet  ini  t  :  » i  med  i  un  sat  a  rat  :  >n ,  '  ru  tendon  ■. 
c  wards  a  linear  dewier  idt  -ru  •••  •  >t’  w  on  t  *.e  pn>|>agat  um  a  is*. am  •••  m-  om<  -  t  . 

‘  n  t*)t  h  oases,  numer  i  ea  .  ,«t  ana  i  >  •  .  •  -a  i  resul's  •  c  >sel  v  agr  -n  v-ii  _  '  n* 

v  a  i  i  d  i  t  v  of  >ir  appr  i\  i  ma  t  .  n  . 


w  ■  ;  ,  a  r  i ,  ,i  an-  *  tier  :  n<)  •  "it 

111,  O''.  e  I  gt  Ti-.  1  .  I*  •  -  arid  •  t  t  •  ■ '  •-  ■ '  V.ll'K  si;  1  t  t  . 

'll  1  r  I  ’  tr'-V'  ases  .  VCal  i  ",  ’  ;  e  ■  1-  les  t  .  -r  I 

of, .si  a,  i  v  reason,  it..-  .  >  ,  i  ,  -  ndt'  ,  "lie  tie's!.-. 
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Fig  1  Temporal  width  of  solitcn  (divided  by  initial  width) 
as  a  function  of  propagation  distance  X,  for  r^-2 
and  y=G.1.  Numerical  integration  (solid)  and 
analytical  results  (dotted). 
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